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Abstract. The present paper establishes a certain duality between the Dirich- 
let and Regularity problems for elliptic operators with /-independent complex 
bounded measurable coefficients (f being the transversal direction to the bound- 
ary). To be precise, we show that the Dirichlet boundary value problem is solv- 
able in V , subject to the square function and non-tangential maximal function 
estimates, if and only if the corresponding Regularity problem is solvable in LP . 
Moreover, the solutions admit layer potential representations. 

In particular, we prove that for any elliptic operator with /-independent real 
(possibly non-symmetric) coefficients there exists a p > 1 such that the Regular- 
ity problem is well-posed in LP, 



1. Introduction 

We consider a divergence form elliptic operator 

L := -divA(x)V, 

defined in R" +1 = {(x, t), x e R n , t > 0). Here A is an (n + 1) x (n + 1) matrix of 
bounded, complex-valued, ^-independent coefficients, which satisfies the uniform 
ellipticity condition 

«+i 

(1.1) m 2 < Re<A(x)^> := Re ^AyC*)^, ||A|| L » (]r) < A~\ 

£.7=1 

for some A > 0, and for all £ 6 C" +1 , x e W. The number A in (1.1) will be 
referred as the ellipticity parameter of L. As usual, the divergence form equation is 

1 2 

interpreted in the weak sense, i.e., we say that Lu = in a domain Q. if u € W^ c (Q) 
and 

(1.2) JJ AVwW = 0, 

for all complex valued ¥ e C^°(fl). For us, Q. will be a Lipschitz graph domain 

(1.3) Q,^ ■- {(x, t) € R" +1 : t > i/r(x)} , 
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where \J/ : R" — > R is a Lipschitz function, or more specifically (but without loss of 
generality), Q will be the half-space R" +1 := {(x, t) e R" x (0, oo)}. We shall return 
to this point below. 

We say that the Dirichlet problem (D p >) for L is well-posed for some 1 < p' < oo 
if given any / e LP (R") there exists a unique solution to the boundary problem 



(*V) 



( Lu = in R% +1 
lim^o u — f n.t. 

JI^v*(")IIlp'(r«) < 00 • 



and 

(1.4) \\NM\\iP' m < Q\f\\LP\w>y 

Here, the notation "lim,_>oM = /n.t." means that lim^-^o) u(y, t) = f(x), for 
a.e. x e R", where the limit runs over (y, t) e T(x) := {(y, t) e : \y - x\ < t}, 
and 

N t F(x)= sup |F(z,0|. 

(z,f)er(.v) 

The Regularity problem (R p ) for L is well-posed for some 1 < p < oo if given any 
/ e Ly(W) there exists a unique solution to the boundary problem 

' Lu = OinRf 1 
(R p ) < lim^oM = /n.t. 

JI^V(Vm)||l/'(R") < °°> 

and 

(1.5) IIAWHiPCR*) < C||V||/|| LP(r) . 

The homogeneous Sobolev space Lj(R n ) is the completion of with respect to 
the Sobolev norm ||V||/||//>(r«). While fairly evident here, it will be convenient to 
distinguish the gradient in and the gradient in R" throughout the paper, and 
we shall denote the latter by Vy .The modified non-tangential maximal function is 
given by 



i 

NF(x) = sup iff \F(y, s)\ 2 dyds 

(z,t)eT(x) \ J J 
\W(z,t) 

where W(x, t) = A(x, t) x (t/2, 3t/2), A(x, f) = {y e W : \x - y\ < t). 

We note, furthermore, that (D p >) and (R p ) above are defined in An analo- 
gous definition applies to the lower half-space. Well-posedness in R" +1 stands for 
the well posedness both in R" +1 and R" +1 , and similarly for other properties. 

It has been proved in [HKMP] that for any elliptic operator L with real bounded 
measurable f-independent coefficients there exists p' < oo such that the Dirichlet 
problem (D p >) is well -posed. The purpose of this paper is to establish a certain 
duality between the Dirichlet and the Regulatity problems, and in particular, to 
show that for any elliptic operator L with real bounded measurable f-independent 
coefficients there exists a p > 1 such that the Regularity problem (R p ) is well- 
posed. Before stating the main result, let us introduce some relevant terminology. 
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Here and throughout the paper, the capital letters X, Y, Z denote points in K" +1 
and the corresponding small ones stand for the points in R". Furthermore, B = 
B R (X) = B(X,R) is the ball in R' ,+1 centered at X e R" +1 with the radius R > 0, 
and A - A R (x) = A(x,R) is the ball in R" centered at x e R" with the radius R > 0. 
Then the tent regions are 7\A) = A = {(x,t) e Wf l : dist(x,A c ) > t}, and the 
Whitney cubes are, as above, W(x, t) = A(x, t) x (t/2, 3t/2), (x, t) <= R^ +1 . 

Throughout the paper L will be an elliptic divergence form elliptic operator with 
bounded, measurable, complex-valued, ^-independent coefficients. We shall as- 
sume, in addition, that the solutions to Lu = in Rl +1 are locally Holder contin- 
uous in the following sense. Assume that Lu - in R£ +1 in the weak sense and 
B 2R (X) c R^. +1 ,X e Rf \ R > 0. Then 



for some constants ji > and C > 0. In particular, one can show that for any p > 



We shall refer to property (1.6) by saying that the solutions (or, slightly abusing 
the terminology, the operator) satisfy the De Giorgi-Nash-Moser bounds. Respec- 
tively, the constants C and fi in (1.6), (1.7) will be referred to as the De Giorgi- 
Nash-Moser constants of L. Finally, following [AAAHK], [HMiMo], [HMM], we 
shall normally refer to the following collection of quantities: the dimension, the 
ellipticity, and the De Giorgi-Nash-Moser constants of L, L* collectively as the 
"standard constants" . 

We note that the De Giorgi-Nash-Moser bounds are not necessarily satisfied for 
all elliptic PDEs with complex ?-independent coefficients [F, MNP, HMMc]. How- 
ever, (1.6), (1.7) always hold when the coefficients of the underlying equation are 
real [DeG, Na, Mo], and the constants depend quantitatively only upon ellipticity 
and dimension (for this result, the matrix A need not be ^-independent). Moreover, 
(1.6) (which implies (1.7)) is stable under small complex perturbations of the co- 
efficients in the L°° norm (see, e.g., [Gi], Chapter VI, or [Al]). Thus, in particular, 
(1.6)-(1.7) hold automatically, e.g., for small complex perturbations of real elliptic 
coefficients. We also note that in the ^-independent setting that we consider here, 
the interior the De Giorgi-Nash-Moser/M estimates hold always when the ambient 
dimension n + 1 = 3 (see [AAAHK, Section 11]). 

The solutions to elliptic operators with real coefficients (and their complex per- 
turbations) also satisfy the boundary Holder continuity condition: 

Assume that Lu - in T(A2r(x)), x e E", R > 0, and u = on A2r(x). Then 



(1-6) 




(1.7) 




(1-8) 




for all (y, s) e T(A R {x)), 
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for some p > and C > depending on the dimension and the ellipticity pa- 
rameter of L only. It is not clear, however, whether the interior Holder continuity 
implies boundary Holder continuity for a general elliptic operator with bounded 
measurable complex coefficients. 

Let us now turn to the layer potentials. Let L be an elliptic operator with bounded 
measurable coefficients. By E, E* we denote the fundamental solutions associated 
with L and L* , respectively, in R" +1 , so that 

L XJ E(x, t; y, s) = 6 (y%s) (x, t) and L* y s E*(y, s; x, t) = L* s E(x, t; y, s) = 6 {xA (y, s), 

where 5( XJ ) denotes the Dirac delta function at the point (x, t). One can refer, e.g., 
to [HK] for their construction and properties. We note for future reference that 
when the coefficients of the underlying matrix are ?-independent, 

(1.9) E(x,t;y,s) = E(x,t-s;y,0), 

and hence, in particular, one can swap the derivatives in t and s for the fundamental 
solution. 

The single layer potential and the double layer potential operators associated 
with L are given, respectively, by 

(1.10) Sff(x) = J E(x,t;y,0)f(y)dy, t e R, x e R", 

R" 

D^f(x) = J d VA ,, y E*(y,0;x,t)f(y) dy, t*0,xe R". 

Here, the conormal derivative is roughly d VA u = -e n+ \A(y)Vu, e n+ i = (0, ...,0, 1). 
The precise meaning of the latter on the boundary will be discussed later, with the 
Preliminaries. 

Going further, it will be convenient to separate the issue of well-posedness of 
(D p '), (R p ) as defined above, from solvability. We shall say that (D p >) (resp., (R p )) 
is solvable if for every / e C™(R n ) there exists some function u € w}£(Wl +1 ) 
such that (Dp') (resp., (R p )) is valid and (1.4) (resp., (1.5)) holds. Throughout 
the manuscript, the constants in (1.4), (1.5) will be referred to as the solvability 
constants of L. We say, furthermore, that (D p >) (resp., (R p )) is well-posed with 
the Green formula representation if for every / £ L p ' (resp., / £ L^) there exists 
a unique u e wj^(R n + +l ) such that (D p >) (resp., (R p )) is valid, (1.4) (resp., (1.5)) 
holds, and 

u = -Df(f) + Sf(d VA u). 

For reasonably nice solutions the formula above is simply the result of inte- 
gration by parts. However, in the context of the general elliptic operators with 
complex coefficients, it seems to be necessary to require its validity explicitly in 
order to establish uniqueness. See Section 7 for a detailed discussion. 

The main result of this paper is as follows. 

Theorem 1.11. Let Lbe a divergence form elliptic operator with bounded, complex- 
valued, t-independent coefficients. Assume, in addition, that the solutions to Lu - 
in Wl +1 satisfy the De Giorgi-Nash-Moser bounds, and that the same is true for the 
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adjoint operator L*. Then there exists s > depending on the standard constants 
only such that for any l<p<2 + s, - + y — I, the following are equivalent: 

(a) The Dirichlet problem (D p >) for L* is solvable in R" +1 , and, inaddition to(lA), 
the solution satisfies the square function bounds 

(1.12) m(^u)\\ I y m <C\\f\\ I y m , 
where 3\ stands for the square function, that is, 

<u3) ™ ixy - = (ILj^P)' 12 ■ 

for F : R n ± +l -> R. 

(b) The Regularity problem (R p )for L is solvable in R+ +1 . 

(c) The Regularity problem (R p )for L is solvable in R+ +1 , and the solution can be 
represented by means of layer potentials, that is, 

(1.14) u(x,t) = Sf (S^y 1 f(x), (x,t) eR" ± +l . 

In particular, the operator f i— > NiVS^f) is bounded in L p (R n ) and 

S£ : Z/(R") -> L((R") 

is invertible. 

(d) The Dirichlet problem (D p >)for L* is solvable in R" , and the solution can be 
represented by means of layer potentials as 

(1.15) u(x,t) = Sf (stf) 1 f{x), ix,t) € Rf 1 . 

In particular, f i— > NiSf f) is a bounded operator from Z/^R") to LP (R n ) 
and 

: z/^R") -> z/(R") 

is invertible. Here L^(R") := (Z,((R n ))*. 

77je solutions to (D /7 <) arcc? (i? p ) in (a) - (J) are unique in the class of solutions 
with the Green formula representation and, respectively, the corresponding bound- 
ary formulas are well-posed with the Green formula representation. 

In combination with the results in [HKMP], [AAAHK], [HMiMo], [HMM], 
and [R] (as regards the latter, see also [GH] for an alternative proof), Theorem 1.11 
yields the following Corollaries. Their proofs can be found in Section 8. 

Corollary 1.16. Let L be a divergence form elliptic operator with bounded mea- 
surable t-independent coefficients, and assume that there exists an operator Lq, 
falling under the scope of Theorem 1.11, such that for some p > 1 one of the 
properties (a) - id) is satisfied for Lq. Assume, furthermore, that 

\\A - Aolto»(R«) < ?/ 

for a sufficiently small n > depending on the standard constants and the involved 
solvability constants ofLo, Lq. Then all four assertions (a) — id) of Theorem 1.11 
are valid for the operator L as well. 
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Corollary 1.17. Let L be a divergence form elliptic operator with bounded mea- 
surable coefficients, and assume that there exists an operator Lq, falling under the 
scope of Theorem 1.11, such that for some p > 1 one of the properties (a) — (d) is 
satisfied for Lq. Denote 



where W{x, t) = A{x, t) X {t/2, 3t/2) and A(x, t) = [y e W : \x - y\ < t}. Suppose 
further that e satisfies the Carleson-type condition 



for a sufficiently small n > depending on the standard constants and the involved 
solvability constants of Lq, Lq. Here Rq := Q X (0,^(2)) is a Carleson box and 
£(Q) is the side-length of the cube Q c R". 

Then all four assertions (a) - (d) of Theorem 1.11 are valid for the operator L 
as well. 

Corollary 1.19. Let L be either a divergence form elliptic operator with real (pos- 
sibly non-symmetric) bounded t-independent coefficients, or a perturbation of such 
an operator, in the sense of Corollary 1.16 or in the sense of Corollary 1.17. Then 
there exists p > 1, depending on the dimension and the ellipticity constant of L 
only, such that statements (a) — (d) of Theorem 1.11 are valid for L. Moreover, the 
corresponding solutions to the Dirichlet and Regularity problems are unique. 

In particular, the Regularity problem is well-posed, for some range of p > 1, for 
any elliptic operator with real t-independent coefficients and for its perturbations. 

We note that this result is sharp, in the sense that one cannot specify the range of 
well-posedness of boundary value problems for real non-symmetric £-independent 
operators which would not depend on the ellipticity parameter of the operator. 
More precisely, for every p > 1 there exists an elliptic operator L with real non- 
symmetric ?-independent coefficients such that the Regularity problem (R p ) is not 
well-posed. A similar statement holds for the Dirichlet problem: given any p' < oo 
there exists an elliptic operator L with real non-symmetric £-independent coeffi- 
cients such that the Dirichlet problem (D p >) is not well-posed. The counterexample 
can be found in [KKPT] and [KR] in the context of the Dirichlet and Regularity 
problem, respectively. 

In the case of real coefficients, the solvability of the Regularity problem for some 
LP, p > 1, is equivalent to the solvability in Hardy spaces, much as the solvability 
of the Dirichlet problem in LP for some p' < oo is equivalent to the solvability 
of the Dirichlet problem in BMO. For the Dirichlet problem the equivalence was 
established in [DKP] and for the Regularity problem in [DK]. We refer the reader 
to [DKP], [DK] for precise statements. Here we just point out that in the realm 
of real coefficients the results of Corollary 1.19 automatically extend to Hardy and 
BMO spaces for Regularity and Dirichlet problems, respectively. 

Finally, we remark that all the results, in particular, Theorem 1.11 and Corol- 
laries 1.16-1.19, automatically extend to Lipschitz domains as defined in (1.3). 



e{x, t) := sup {\A(y, s) - A (y)\ : (y, s) e W(x, t)} 



(1.18) 
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This is a consequence of the fact that a "flattening" change of variables (x, t) h> 
(x, t - ijj{x)), which maps a Lipschitz domain {(x, t) e : ? > tA(x)} into 
preserves the class of f-independent elliptic operators. 

Let us now discuss the history of the problem. The study of elliptic boundary 
problems (D p >), (R p ) has started with the results for the Dirichlet problem for the 
Laplacian on Lipschitz domains [Dl], [D2]. The first breakthrough in the context 
of the elliptic operators with bounded measurable coefficients came in [JK], where 
the authors realized how to resolve (D p >), p' - 2, resting on the so-called Rellich 
identity. The latter is essentially a result of an integration by parts argument which 
allows one to compare the tangential and normal derivatives of the solution on the 
boundary, in the sense that 

(1-20) II v ii/IIl 2 (r«) ~ II^v a «IIl 2 (R")> / = u \ R n- 

The Rellich identity underpinned the development of the elliptic theory for real 
and symmetric operators, and over the years the problems (D p >) and (R p ) were 
resolved for the sharp range of p for operators with real symmetric £-independent 
coefficients in [KP] and perturbation results in the spirit of Corollary 1.17 were 
obtained in [D3], [FJK], [FKP], [KP2], (see also [AAAHK], [AA], [HMM] for 
later developments in connection with the perturbation questions). 

We remark that some "smoothness" of the underlying matrix in t is necessary 
for well-posedness [CFK] and thus starting the investigation with the £-independent 
case is natural in this context. 

The argument for the Rellich identity heavily used the condition of the symme- 
try of the matrix, and thus, could not be extended neither to real non-symmetric, 
nor more generally, to the complex case. The only exception to this rule was the 
resolution of the Kato problem [CMcM, HMc, HLMc, AHLMcT], in which (1.20) 
was established in the absence of self-adjointness, in the special case that the ma- 
trix has block structure, that is, A = {A^}"^ with Ay )JJ+ i - A n+ ij = 0, j = 1, n. 
The observation that the solution of the Kato problem amounts to (1.20) is due to 
C. Kenig, see [K]. Moreover, quite recently (simultaneously with the preparation 
of this manuscript) it was shown in [AMM], by a refinement of the proof of the 
Kato conjecture, that the Kato estimate, or more precisely, the > side of (1.20), 
could be extended to the "block triangular" case when only Aj tn+ \ = without 
necessarily A n+ \j = and similarly, the < side of (1.20) holds when A n+ ij = 0. 

However, the aforementioned ideas could not be directly applied to the general 
case of a non-symmetric matrix lacking any additional block structure 1 . Moreover, 
it was demonstrated in [KKPT] (see also [KR]) that the well-posedness in L 2 may 
fail when matrix has no symmetry and thus, (1.20) is not to be expected. Nonethe- 
less, using a completely different approach, in [KKPT], [HKMP] the authors have 
established that for any operator L with real non-symmetric coefficients there is a 
p' < oo such that (D p >) is solvable. This raised the question of solvability of the 
Regularity problem. 



although the technology of the Kato problem continues to play a crucial role in the present paper 
and in [HKMP]. 
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In the particular case of the operators with real non-symmetric coefficients in di- 
mension two and their perturbations the Regularity problem was resolved in [KR], 
[B]. The present paper resolves this problem in arbitrary dimension. It shows that 
the solvability of the Dirichlet problem is generally equivalent to that of the Reg- 
ularity problem and thus, in the context of real non-symmetric matrices there is 
always a p such that (R p ) is solvable. The core of our argument is a new Rellich- 
type inequality. We demonstrate that, in fact, for any operator L with complex and 
^-independent coefficients the solvability of the Dirichlet problem (D p >), together 
with the square function bounds, entails a one-sided Rellich inequality, 

(1-21) \\d VA u\\LP(R») £ I|V||/||ZP(R»), / = u\ Rn . 

This ultimately paves the way to (R p ). Of course, having a reverse inequality as 
well would be extremely interesting, but at the moment seems quite challenging. 

Finally, we also point out that the actual question of connections between the 
Dirichlet and Regularity problem has received considerable attention in the litera- 
ture, and some partial results were established in [V] ((R p ) <==> (D p >), Laplacian 
on a Lipschitz domain), [KP] ((R p ) ==> (D p >), real coefficients), [S], [KiS] ((R p ) 
<=> (D p '), real symmetric constant coefficient systems), [KR] ((D p >),=$» (R p ), 
real coefficients, dimension two). Some related counterexamples were obtained in 
[M]. Our main result, Theorem 1.11, generalizes all implications above, at least 
as far as the ?-independent matrices are concerned. The proof of (1.21) builds on 
Verchota's duality argument [V] , reducing matters to proving LP estimates for 
certain conjugates, and in turn, the estimates for the conjugates will be obtained by 
an extension of an argument in [AAAHK] which exploits the solution of the Kato 
problem. 



2. Preliminaries 

Let L be a divergence-form elliptic operator with f-independent bounded mea- 
surable coefficients. Any solution to Lu = satisfies the interior Caccioppoli 
inequality: 

Assume that Lu = in R'^ +1 in the weak sense and B 2R (X) ctf'je 
R > 0. Then 

(2.1) jj |V M (F)| 2 dY < ^ jj \u(Y)\ 2 dY, 

Br(X) B 2 r{X) 

for some C > depending on the dimension and the ellipticity parameter of L only. 
An analogous statement holds in 

Recall that we assume, in addition, that solutions of L and L* satisfy the interior 
Holder continuity conditions, that is, the De Giorgi-Nash-Moser estimates (1.6), 
(1.7). 

We proceed to the issues of the non-tangential convergence and uniqueness of 
solutions. Unless explicitly stated otherwise, we assume throughout the rest of the 
paper that L = - div(AV) is an elliptic operators with complex bounded measurable 
£-independent coefficients and that the solutions to Lu = and L*u = satisfy the 
De Giorgi-Nash-Moser bounds. Furthermore, throughout the paper we assume that 
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n > 2, as much of this theory in the case n - 1 has already been treated in [KR] 
and [B]. 

Lemma 2.2. Suppose that u e W^(R" +1 ) is a weak solution of Lu = 0, which 
satisfies N(Vu) € L p (R n )for some 1 < p < oo. Then 

(i) there exists f e L^(R") such that u —> f rut. a.e., with 

\u(y, t) - /(x)| < t N(Vu)(x), for every (y, t) e T(x), x e R n , 
and 

\\f\\ LI ;<\\N(Vu)\\ LP ; 

(ii) for the limiting function f from (i), one has 

V|| M (-,0 — > V|,/ 

as t — > 0, in the weak sense in L p ; 

(iii) there exists g e L p (R n ) such that g = 5 Va m in f/je variational sense, i.e., 



ff 



A(X)V«(X)VO(X) - / ^o(x) Jx, 



/or O 6 Cg°(R" +1 ) W := <D| f=0 ; 
(iv) for the limiting function gfrom (iii), one has 

-e n+ iAVu(; t) — > g 

as f — > 0, in the weak sense in LP. Here e n+ \ = (0, 0, 1). Finally, 

(2.3) \\g\\ LP (R«)<C\\N(Vu)\\ LP(m . 

An analogous statement holds in K." +1 . 

The Lemma can be found in [HMiMo] (as stated above), [HMM] (for somewhat 
more general operators), [AAAHK] (for p = 2), [KS] (for real coefficients) and the 
proof in all cases closely follows an analogous argument in [KP]. 

Remark. We note that the convergence results above entail, in particular, the fol- 
lowing. The solvability of (R p ) as defined in the introduction, that is, for C^°(R n ) 
data, implies existence of a solution to (R p ) for any / 6 Lj(R"). An analogous 
statement is valid for the Dirichlet problem. We refer the reader, e.g., to [KS], The- 
orems 4.6, 5.6, where these results were established for real symmetric matrices, 
and the proofs apply to our case without changes. 

The next order of business is to define a weak solution to the boundary problem 
Lu = in R n + +l , u - f € C~(R"). In principle, taking a harmonic extension 

of / to R" +1 (denote it by w) and then using the Lax-Milgram lemma to resolve 
Lu = -Lw with zero trace on the boundary, we get a solution in W 1,2 (R" +1 ), the 
factor space of functions modulo constants with the seminorm given by the norm 
of the gradient in L 2 (R" +1 ). It is somewhat more convenient though to use a non- 
homogeneous space. One option (and here we follow an approach in [KR]) is to 
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work in the following framework. Let W l,2 (Wl +l ) denote the space of functions F 
for which 

One can define the trace operator, for instance, as Tr : W 1,2 (R" +1 ) — > L 2 (R n ), 
a continuous extension of the restriction to the boundary operator, with L 2 (R") 
denoting the space of functions / on R" with 

ii/ifc«, = (/ E ,i/u)i 2 T ^)" 2 <»- 

(follow, e.g., the argument in [E], p. 272). Notation W ' (R" +1 ) stands for the space 
of functions in W l ' 2 (B," + +l ) with trace zero. 

Lemma 2.4. For any 2 < r < 2( £*^ and I < q < n > 1, and for any 

f e L 2 (R") n L'\R n ) n L q (W) there exists a unique u £ W U2 (R n + +i ) such that 
Tr u - f and Lu - 0, in the usual weak sense. Moreover, 

(2.5) IMIjyU^-H, £ ll/llz,2(R») + II/IIl'(R") + II/IIl«(R")- 

Proof. The proof is a modification of an analogous arguments [KR]. Here we only 
mention the main idea. The remaining details are quite easy to fill in, and if needed, 
the reader may consult [KR]. 

The basic idea, already mentioned above, is to realize u as v + w, where w is the 
solution to the Laplace's equation with data /, that is, the Poisson extension of /, 
and v is the solution to Lv = G (where G = -Lw) with zero boundary data, given 
by the Lax-Milgram Lemma. A direct computation shows that w e W i ' 2 (E" + +l ) 
satisfies (2.5). The Lax-Milgram lemma assures existence of the unique solution 
v e Wo'\M" + +1 ) to the problem 



(f AVvV®dxdt= II AVwVOdxdt, for all O € W U (I 



r 1 ). 



This requires boundedness and coercivity of the bilinear form with respect to the 
norm in IV'' 2 (R" +1 ). Boundedness is obvious from the definition, and the coercivity 
follows from the Poincare inequality. We note for the future reference that the Lax- 
Milgram lemma, in particular, assures that 

(2.6) IIVvll^^+i) < C\\Vw\\ L 2 (K+ i y 

Finally, the desired estimates (2.5) follow from the combination of aforemen- 
tioned bounds on w and (2.6). □ 

Let us now provide somewhat more precise asymptotic estimates on the weak 
solution u constructed above in the case when / is a nice compactly supported 
function. 

Lemma 2.7. For any f e Cq(W) the weak solution u e W 1>2 (R" +1 ) to the problem 
Lu = 0, u = f, warranted by Lemma 2.4, satisfies 

(2.8) \u(x, 01 < C f t 1 "^, (x, t) G R n + + \ 
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provided that t is significantly bigger than the size of the support of f. The constant 
C depends on f and the operator L. 

The statement and the proof of this Lemma is the only place in the paper where 
constants denoted by C are allowed to depend on data /. The results, however, will 
be only used qualitatively to ensure the convergence of some later arising integrals. 

Proof. Recall the construction of the solution from the proof of Lemma 2.4. The 
classical harmonic Poisson extension of /, denoted by w, clearly satisfies (2.8). 
In fact, it decays faster, as w(X) = 0(\X\~ n ), at infinity. It remains to estimate 
v. To this end, we decompose v as a sum £°V;, where v ' e Wq'\WI +1 ) is the 
Lax-Milgram solution to 



1,2/qnin+l-j 



\\ AVvjVQdxdt = \\ AVwiVQdxdt, for all Ueff^ 

Here w ; - = rjjW, with ijt, i = 0,1,..., being the elements of the usual partition 
of identity associated to dyadic annuli of radius 2' centered at the origin (and 770 
associated to the unit ball). Then, in particular, for all i sufficiently large depending 
on the size of the support of / we have 

(2.9) l|Vv,-|b (Rrl) < CllVwdb^i) < C2-''(" +1 > /2 , 

using asymptotic estimates for the harmonic functions. 

Next, observing that v, = on the boundary we employ Poincare inequality to 
write 

v 1/2 / v 1/2 

ff\vi\ 2 dx\ <C2 k I jj\Vv£dx\ < C 2K^)2-^I\ 

\B k (0) I \B k (0) I 

Finally, by De Giorgi-Nash-Moser estimates, 

/ v 1/2 / \ 1/2 / \ 1/2 

\u{x,t)\< jj \u\ 2 dx\ < jj \w\ 2 dX\ +J2 jj hfdX 

\B t /2(x,t) J \B t / 2 (x,t) J ,= 1 \B tn (x,t) 

-n |_5±1 
< t + t 2 , 

with the implicit constant depending on / and assuming that t is large enough 
compared to the size of the support of /. □ 

As discussed above, any u e W l ' 2 (Wl +l ) has a trace in L 2 (R"). Moreover, for 
any u which is a solution to Lu = in the sense of Lemma 2.4, one can define the 
conormal derivative d v u = d VA u, in the sense of distributions, via 



(2.10) dyii(x) <p(x) dx := A(X)V M (Z)VO(Z) dX, 



for any tp e Cq(W) and O e C ( ~(R" +1 ) such that <p := O | f=0 . The details are as 
follows. 
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Lemma 2.11. Suppose that u € W i,2 (M.'l +i ), and that Lu - in R" . Then d VA u 
exists in L 2 _ l/2 (M. n ); i.e., there is a g u e L 2 _ 1/2 (R") such that for all H € W l2 (M. n + +l ), 
with trace tr(H) = h € L 2 / 2 (R"), we have 



where {•, •) denotes here the duality pairing of L_j , 2 (R") ow<i ^-mO^")- The analo- 
gous statements hold for the adjoint L* , and in the lower half-space. 



Proof. We define a bounded linear functional A„ on L^ 2 (IL") as follows. For h € 

L? /2 (R"), set 



where H is any W ' (R" + ) extension of h (of course, such extensions exist by 
standard extension/trace theory). Note that A„ is well-defined: indeed, if Hi and H2 
are two different W ia (Rl +1 ) extensions of the same h, then H l -H 2 € W^ 2 (R'l +i ), 
whence it follows that 



since Lu = in the weak sense. Moreover, it is obvious that A„ is linear. To 
see that the functional is bounded, we simply choose an extension H (e.g., the 
harmonic extension), for which 

H#llw"- 2 (Rr') - C ° IWIi? /2 (K.") ' 
for some purely dimensional constant Co- We then have 

\A u (h)\ < \\A\\ 

i.e., ||A U || < ||w||^i,2 The conclusion of Lemma 2.1 1 now follows by the Riesz 
Representation Theorem. □ 

Lemma 2. 1 1 allows us to justify the definition of the conormal derivative by 
(2.10) for all functions u e W i,2 (Wl +1 ), by first identifying such a function with 
the corresponding equivalence class in W /I ' 2 (R." +1 ), then applying Lemma 2. 11 and 
then deducing that for any particular representative of this equivalence class, in 
particular, for u itself, we have (2.10). Such a definition gives identical result for 
u in the same equivalence class, but, naturally, the conormal derivative would not 
distinguish functions that differ by a constant. 



Note that if, in addition, N(Vu) e L p (R n ), for some 1 < p < 00, then d v u € 
LP(R n ) by Lemma 2.2. 



Theorem 3.1. Let L be an elliptic operator with t-independent coefficients such 
that the solutions to Lu = and L*u = in R+ +1 satisfy the De Giorgi-Nash- 
Moser estimates. Let u be a solution to the Dirichlet problem Lu = in R" , 



(2.12) 






3. Boundary estimate: a version of the Rellich-type inequality 



u 



0F1 



- f, for some f € C^°(R"), in the sense of Lemma 2.4. Suppose that for 



REGULARITY PROBLEM 



13 



some 1 < p' < co the Dirichlet problem (D p >)for the operator L* is solvable with 
the square function bounds, that is, the solution u satisfies both (1.4) and the square 
function estimate (1.12). Then the variational derivative ofu defined by (2.10) can 
be identified with an LP function, and 

(3-2) \\dy A U\\ L „(Rn) < C \\V\\f\\lJ>(«."), 

where - + y — 1 . The constant C depends on the standard constants and on the 
solvability constants ofL* involved in (1.4) and (2.10). 



Proof. We aim to show that for u, a solution to Lu - in R" +1 , the normal deriva- 
tive on the boundary is controlled by the tangential derivatives for some 1 < p < oo. 
To this end, take 1 < p' < oo such that (D p >) for the operator L* is solvable with 
the square function bounds, consider any g e with \\g\\jy < 1, and denote by w 
a solution to (D p ') for the operator L* with boundary data g. Then 

(3.3) / d VA ugdx= \ \ VF(X) A*(X)Vw(X) dX = / fd~wdx, 



= f. This follows simply from the definition 

R" 



for any F € C£°(R" +1 ) such that F 

of the weak conormal derivative in (2. 10). Indeed, by definition 

J fdl~wdx = JT VF(X) A*(X)Vw(X) dX - Jj Vw(X) A*(X)Vw(X) dX, 

where for the second equality we use the fact that u - F € Wq' 2 (R" +1 ) and w is a 
solution in the sense of (1.2) (evidently, functions are dense in W 1,2 (R£ +1 )). A 
similar argument applies to show that 

J d VA ugdx = jj A(X)Vu(X) Vw(X) dX. 

We shall take an extension of / in the form F(x, t) :- f{x)n r n{t), R » r, where 
T] r>R (t) € C^((-(R + r),R + r)), T] r<R (t) — 1 for / G (~R,R), and \t]'\ < l/r. Here r is 
chosen so that A r contains supp(/). Then the right-hand side of (3.3) is equal to 

(3.4) Jj~ VF(X) A*(X)Vw(X) dX + JT VF(X)A*(X)Vw(X)dX 

A r x(0,R) A r X(R,r+R) 

= d jf(x) e j A*(X) Vw(X) dX + Jj VF(X)A*(X)Vw(X)dX 



7=1 



A,x(0,R) A r x(R,r+R) 



= it,f d if^ (J ejA*(x)Vw(x,t)dt^j dx+ VF(X)A*(X)Vw(X)dX, 

j =l A,- ' A r x(R,r+R) 

where ej = (0, 0, 1, 0, 0) is the j-th basis vector of R" +1 , and we used the fact 
that by construction F = f in A r x (0, R) and hence, in this range, it is independent 
of t. 
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We remark that, intuitively, one can think of the functions in parentheses above 
as an analogue of harmonic conjugates. That is, given a solution w to Lw = in 
R" +1 , a system of L-harmonic conjugates could be defined as follows: 

/CO 
ejA(x)Vw(x,s)ds, (x,t) elf , ; = l,...,n+ 1, 

(see, e.g., [FS] and [KR] for analogous constructions in the case of harmonic func- 
tions and variable-coefficient operators in dimension 2, respectively). Thus, (3.3)- 
(3.4) and forthcoming calculations are actually manipulations with harmonic con- 
jugates in disguise. However, because of the weak nature of the available definition 
of solution and conormal derivative, we have to carefully keep track of the error 
terms. 

Going further, departing from the right-hand side of (3.4), we can write 



(3.6) 




First of all, we claim that the terms # and E^^r both vanish as /? — » oo. Indeed, 




using the Cauchy-Schwarz and Caccioppoli inequalities for the first bound and 
Lemma 2.7 for the second one. Analogously, using Lemma 2.7, we see that E\r < 
C/ ;r /? 1_ T" as well. 

It remains to analyze Ir and //. The integral in // directly gives the desired 
bound by ||V||/||iP. The estimate on Ir is trickier. Using, as before, the decay of w 
assured by Lemma 2.7, we see that it is enough to bound 

(3.7) /:= ( V||/(x)-A*(x)V||v(x)^x, 

JR" 

where formally 

/CO 
w(x, s)ds, v(x) := v(x, 0), 

as the error vanishes as R — > oo. Here and in the sequel, Ay denotes an n x n block 
of the matrix A, that is, {Ajk}" k=l , and L\\ = - divy A||V||, interpreted, as usually, in 
the weak sense. 
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Let us discuss the definition of v. First of all, v(x, t) itself is well-defined for any 
t > as an absolutely convergent integral (using Lemma 2.7), and Vv(-, t) belongs 
to Lf oc (R") for any t > 0. Further, Vyv on the boundary is well-defined as an 
Lf 0C (R n ) function (using the fact that w e W h2 (Wl +1 ) and Lemma 2.7), and V||v(-, t) 
converges to V||v(-) in the sense of distributions and in L| )C (R"). This clarifies the 
sense of (3.7)-(3.8). We note for the future reference that d t v = w belongs to 
^i,2( R «+i) and sa ti s fies Lemma 2.7, (2.1), and (1.8). 

We claim that 



(3.9) 



V||/(x)-A*(x)V||v(x)^x 



iiV||/ib [mimv)^ + iia^v)u l ,/) , 



where 9K, as before, stands for the square function (1.13). Estimate (3.9) is one 
of the core components of our approach to the regularity problem and we will 
concentrate on its proof in the next section. It is interesting to point out that the 
particular quadratic form appearing here is important. Indeed, the case of interest 
is p < 2. The function v itself is a solution to the Dirichlet problem for L* , and it is 
not in general true that the derivatives of the solution belong to some LP , p' > 2. 

For now, let us finish the proof of the Theorem assuming (3.9). 

Recall that d t v = w and by definition w is the solution to (D p > ) for L* satisfying 
both the non-tangential maximal function and the square function estimates. All in 
all, then (3.6)-(3.9) guarantee that for any g e LP 

(3.10) ( dytUgdxZWwfWvMrf. 

JR" 

Hence, 

(3-11) \\d VA u\\LP(R,") < l|V||/||iP(iR«), 

with / - u\n>n. This finished the proof of Theorem 3.1, modulo (3.9). □ 



4. Proof of the mapn estimate (3.9) 

In this section, we establish the "main estimate" (3.9) (re-stated as Theorem 4.10 
below) thereby completing the proof of Theorem 3.1. We shall adapt the proof of 
[AAAHK, estimate (5.3)], which is essentially the case p — 2 of (3.9), and which 
exploits the solution of the Kato problem. In our case, we require LP versions of 
the Kato estimate (cf. (4.2) below). 

Let us start by recalling a few results regarding the square roots of elliptic opera- 
tors and the corresponding square functions that will be used throughout the proof. 
Retain the definitions of N*, N and Jl from Sections 1 and 3, and let 



dzdt 



n+l 
-z\<at 
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where, as before, A = {(x, t) : dist(x, A f ) > t] and A(x, t) - {y e R" : \x - y\ < t}. 
The aperture index a will usually be omitted unless it plays an explicit role in the 
proof. Also, as per (1.13), = ffix, and C = C2. 

For a Lebesgue measurable set E, we let M(E) denote the collection of measur- 
able functions on E. For < p, q < 00 we define the following tent spaces: 



(the notation "sup" is interpreted as the essential supremum) and W(x, t) = A(x, t) x 
(t/2, 3?/2). The spaces T|'(R" +1 ), < p,q < 00 were first introduced by Coifman, 
Meyer and Stein in [CMS]. The spaces f£(Ii" +1 ) and T~(R'| +1 ) started appearing 
in the literature more recently, naturally arising for elliptic PDEs with non-smooth 
coefficients. 

As usual, we say that a family of operators {T t } t >o satisfies L p - L q off-diagonal 
estimates, 1 < p < q < 00, if for arbitrary closed sets E,F cl" 



for every t > and every / € L p (M, n ) supported in E. We remark that whenever L 
and L* both satisfy the De Giorgi-Nash-Moser property, the heat semigroup f t := 
e -f2L|1 , t > 0, satisfies pointwise Gaussian upper bounds, and hence, L p - L q off- 
diagonal estimates for all 1 < p,q < 00. Indeed, if the solutions of L have the De 
Giorgi-Nash-Moser bounds, then so do the solutions to Ly. To see that, let A = 
A(x, r) be an ^-dimensional ball, and let B = B((x, r), r/2) be the corresponding 
(n + l)-dimensional ball. Let u - u{x) solve L\\u = in 2A. Set U(x, t) := u(x), so 
that U is ^-independent. Since U and also the coefficients of L are ^-independent, 
we have that LU(x,t) = L\\u(x) = for all (x,t) e 2B. Since U satisfies the 
De Giorgi-Nash-Moser property in B, we conclude that u satisfies the De Giorgi- 
Nash-Moser property in A. Furthermore, if the solutions to both L\\ and Ljj have 
the De Giorgi-Nash-Moser bounds, then the heat kernel, that is, the kernel of the 
semigroup P t :- e~ n, t > 0, satisfies pointwise Gaussian upper bounds, and 
enjoys Nash type local Holder continuity, by [AT] (Theorem 10, p. 34, loc. cit.), 
as desired. 





where in the last definition, for F € M(R" +1 ), we set 

s(F)(x,t):= sup \F(y,s)\ 



(y,s)eW(x,t) 



(4.1) 
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The latter, in turn, imply that the square root estimate, 



(4.2) 



hf 



<C||V/|| iP , 



holds for all 1 < p < oo, / e Cg°, (the case p = 2 corresponds to the Kato problem 
solved in [AHLMcT], and the generalization to other values of p (given the L 2 
case), in the presence of Gaussian bounds, can be found in [AT]). The Gaussian 
bounds also imply that 

(4.3) \me t f)\\ LP < qi/ib, 

with 9, = t^L\<P t ,Q t = f 2 V|, JQP U 0, = t 3 Ll /2 P„ 6, = f 4 V,|Lj /2 P t , 

holds for all 1 < p < oo, / e LP . In the required generality these square function 
estimates do not seem to be explicitly stated anywhere, but they are all essen- 
tially well known. Indeed, one may verify the case p = 2 by a standard "quasi- 
orthogonality" argument; for p > 2, one may follow the well known argument of 
[FS] to get a Carleson measure estimate when / e L°°, and then use tent space in- 
terpolation to obtain all p e [2, oo); for p < 2, one may first prove the Hardy space 
bound ||>?l(0f/)||£i < CH/lljyi, by a standard argument using the atomic decomposi- 
tion of H { , and the local Holder continuity of the heat kernel, and then interpolate 
to get the full range of p. We omit the details. 

Going further, we record the following result essentially following from the 
Poincare inequality. It was proved for p = 2 in [AAAHK] (Lemma 3.5, loc. cit). 



Lemma 4.4. Assume that a family of operators {R t \t>o satisfies Lr—Lr off-diagonal 
estimates and that R t \ =0 (in the sense o/L 2 oc (R")j. Then 

(4.5) \w£tF)\\» < c\\mw)\\^ 

for every F with tV\\F e T^, i.e., such that the right-hand side of (4.5) is finite, and 
for every 1 < p < oo. 

In fact, a weaker off-diagonal decay rate of R t than the exponential estimates 
above would suffice, but for all relevant choices of /? f 's in the present paper the 
exponential decay will be valid. 

Proof Let us denote A,,, = {y € R" : \x-y\< t], S j(A Xtt ) = A x2 i+i t \ A x , 2Jt , j e N, 
and F x t := f. F(y, t) dy. Then 



, 1/p 

1, 



(4.6) m l (R t F)\\u, = ( / ( || JR tF (y,t)\ 2 ^pj P/ dx 





\Rt([F(; t) - F xat ]XA x2t )(y)\ 2 ^f) Pl dx 
\x-y\<t 1 J j 



Up 



I! 

JJ\x-y\<t 



\R t ([F(; t) - F x at]xsj(^))(y)\ 2 ^P\ " dx 



Up 



I + 11. 
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Using the uniform in t boundedness of R t in L 2 (R") (following from the L 2 - L 2 
off-diagonal estimates) and then the Poincare inequality, we deduce that 



On the other hand, by L 2 - L 2 off-diagonal estimates, 
(4.8) 



7=1 



R" \ J 



e 



r a \ p/1 \ UP 

-(2W) / lF(yj) _ Fx2tl 2 dy AL\ dx \ 

JSj(A x , t ) 1 J J 



where we used the representation 



F(y, t) - F x2t = F(y, t) - F Xt2 j + i t + F x<7j+ i t - F xXn + ... + F x 2 i t - F x2t 
for the second inequality. Using now the Poincare inequality, we have 
(4.9) 




where M > n/2 can be arbitrarily large constant. However, (4.9) simply says that 

CO 

n<Y. 2 ~ Mk 

which in turn implies that 

CXI 

i+u < Y,^ MUC{n ' p)k \\^\^\\F)\\D' < ii^cfViiioiiip, 

k=0 

since M can always be taken large enough. We remark that the sharp constant 
C(n, p) appearing in the change-of-aperture square function estimates was obtained 
in [A2], although for the purposes of the present argument we only need to know 



M 2k+ \tV ]{ F) 
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that the dependence on the aperture is polynomial, and this was already established 
in [CMS]. □ 

At this point we are ready to turn to the proof of estimate (3.9). 

Theorem 4.10. Assume that L is an elliptic operator with t-independent coeffi- 
cients, and that L and L* satisfy the De Giorgi-Nash-Moser bounds. Let w be the 
weak solution to the Dirichlet problem for L* with some data, guaranteed by 
Lemma 2.4, and v be its antiderivative, defined by (3.8). Then 



(4.11) 



V ll f(x)-AUx)V ll v(x,0)dx 



< cnvn/ib> (umWtvMv,' + m*(d t v)\\iy) 



for every f e C^K.' 1 ) and 1 < p < oo. Here C > depends on the standard 
constants only. 

Clearly, (4.11) is only of interest when the right-hand side of (4.11) is finite. It 
is the case, e.g., when the Dirichlet problem is solvable in LP , and the solution 
satisfies the square function estimates. However, this information is not needed to 
establish (4.11). 

Remark. For future reference, we point out that the argument will establish a more 
general result. To be specific, for L and / as in the statement of Theorem 4.10 we 
shall demonstrate that 



(4.12) 



V||/(x)-A*(x)V||v(x,0)dx 

< C||V||/|b> (\mtWtv))\\ LI >' + WNMvJWu,' 



+ 



ff 



dfP t f(x)L*v(x, t) tdtdx 



for every v 



ii/l+l 



with reasonable decay properties sufficient to justify con- 



vergence of involved integrals. Here P t 



, t > 0, is, as before, the heat 



semigroup associated to the operator Ly. Clearly, when v is a solution, as in the 
statement of Theorem 4.10, the last integral on the right hand side of (4.12) is 
equal to zero and (4. 12) reduces to (4. 11). 

Proof of Theorem 4.10. To begin, 



(4.13) 



V||/(x)-A*(x)V||v(x,0)Jx 



-K 



„+i It ( V H^/^ -A*(x)V||v(^,0) dtdx 



= -ff 1 V||3r 
JJr" +1 



P,f(x)-AT,(x)V n v(x, t) dtdx 



JJr" +1 



f(x)-A?,(x)¥\\d t v(x, t) dtdx 



ff 



V\\P,f(x) ■ A*(x)V||^v(x, t) tdtdx 
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+ 2 II V\\d t Ptf(.x) ■ A|(jc)V||5,v(jc, t)tdtdx 



+ 



ff 

ff ,V||^ t . 



f(x) ■ A?,(x)V\\v(x, t) tdtdx =:I + II + III. 



The convergence of integrals and integration by parts above and through the argu- 
ment is justified by our assumptions on v (in particular, the properties of v and d t v 
outlined in the paragraph above (3.9)), and off-diagonal estimates for f t and its 
derivatives. 



Then 
(4.14) |/| = 



II 



L\\'Ptf(x)dfv(x, t) tdtdx 



M(t ^L\{P t { y/hfrf)) 



LP 



Imtfv)^ < llVy/lb \\M(td*v) 



using boundedness of the square function based on t \/L\\P t (4.3) and then the Kato 
estimate (4.2), both available for 1 < p < oo. Largely by the same argument, 



(4.15) |//|< 



^l(f 2 V|| sfQPti y/L^f)) H^Vii^v)!^, 



LP 



< 



l|V||/|b || m^\\d t v)\\^ . 



The estimate on III is more delicate. We write 
(4.16) 111 



= ff f t P 



f(x) L*v(x, t) tdtdx 



+ J2 ff d^ t f(x)djA* n+l d t v(x,t) tdtdx 
n+1 „„ 

+ V I I d}p t f(.x)A* +1 IfjdMxJ) tdtdx -. Ilh + II h + Hh- 
~[ J J '•' 

The term III\ is left alone for the moment. It is zero when v is a solution and 
otherwise, it shows up explicitly on the right hand side of (4.12). The term IIIj 
can be handled just as I and II above and gives the bound 

(4.17) in, < HVy/ib \mm t v)\\^ , 

once again using (4.3) and (4.2). 
As for III 2 , 

(4.18) /// 2 = CV ff djd 2 t <P t f(x)A* +l (I - P t )d t v(x, t) tdtdx 

n _ _ 

+ C Y1 I I , ^tfix) (djAl n+1 P t d t v(x, o) tdtdx =: 7// 2>1 + 7// 2 , 2 , 
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where P t is a nice approximation of identity, e.g., the heat semigroup of the Lapla- 
cian. Since R t :- I — P t kills constants, and, by virtue of standard heat kernel 
bounds, satisfies off-diagonal decay estimates, one can use the Poincare inequality 
to produce the gradient and get an estimate on llh,\ akin to that for I + II. Indeed, 
according to Lemma 4.4, square function estimates (4.3) and Kato estimate (4.2), 



Concerning 77/2,2 , we further write 
(4.19) 777 2 , 2 



=±f{ 



4t 2 L 2 e~ t2L " - 2L\\e- t2 ^ J f{x) (djA* jn+l P t d t v(x,t)j tdtdx 



C^L\e-^ +C 2 L ll e-^ L A f(x)P* t/2 (djA* n+l P t d t v{x,t)) tdtdx 



tff 



C^Lje'^ 2 ^ + C 2 L\\e~* ,2L ^ f(x)x 



x [Vj A ln + l p tdtv(x,t) - (p; /2 djA*. n+i )(P t d t v(x,t))) tdtdx 

n „„ 

JJ R „ +1 (c^Lje-^ + C 2 L l] e'^ L ^ f(x)x 



x (Pt/2 d J A ln + t> ( p td t v(x, 0) tdtdx = U1' 22 + 7//£ 2 , 

where (P*/ 2 djA*j n+l ) (P t d t v(x, t)) is interpreted as a product of two functions, while 
< P*i2^A*jn^\ p tdtv{.x, t), as before, is a result of an operator 'P*/ 2 djA* n+l P t acting on 
the function d t v. 

Much as in the analysis of III%\, we note that the operator 

R t := t {P; /2 djAl n+1 P t - (P; /2 djAl n+1 )P t ) 

kills constants and satisfies L 2 - L 2 off-diagonal estimates (see, e.g., [AHLMcT]), 
so that by Lemma 4.4 one can bound 777 2 2 similarly to Finally, 

(4.20) III'i 2 < II^^AWIL x 

x |^ ((Ci? 3 LV?' 2L + C 2 tLe-^' 2L ^ f) || \\N*(Ptd,v)\y . 

The latter estimate follows from the duality of tent spaces T 2 and T 2 and the fact 

that r 2 °° • Tt ^ T 2 (see [CMS] for the case p' > 2 and [HMM] or [CV], p. 313 
for all < p' < oo). 
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For each j, the estimate on the Carleson measure of C{tP*/ 2 djA*j n+l ) is con- 
trolled, since the Kato problem for L| is equivalent to the statement that 

||C(#>; /2 V|| • £)lk-(R») < C||£||l» ( r»), 
for every t> e L°°(R'\ C") [AHLMcT]. Hence, 
(4-21) < ||V||/||„ \\N*(d,v)\y ■ 

Indeed, recall that P t is a "nice" approximation of identity, e.g., the heat semigroup 
of the Laplacian and hence, one can drop it inside the non-tangential maximal 
function using the fact that N t .(P t (F(-, i)) < M(N*F) with M denoting the Hardy- 
Littlewood maximal function. □ 

This finishes the proof of Theorem 3.1. 

5. BOUNDEDNESS OF LAYER POTENTIALS AND ESTIMATES FOR SOLUTIONS IN 

Let us start by recalling a few known estimates on the layer potentials and related 
operators that will be used throughout this section. 

Proposition 5.1. Suppose that L is t-independent, that L and its adjoint L* satisfy 
the De Giorgi-Nash-Moser bound (1.6). Then the single layer potential satisfies 
the following estimates in L 2 : 

, L i2 dxdt 2 



(5-2) / / td t VSff(x)\ — < CU/. 

J -co JB." \t\ 

(5.3) sup (|| (S? V) /|| L 2 (ro + ||V5f/|| L 2 (Rn) ) < \\f\\ L 2 m , 



(5-4) \\N(VS L ±t f\\ L 2 (Rn) <\\f\\ L 2 (Rn) . 
Moreover, there exists e > such that for all 1 < p < 1 + s we have 
(5.5) l|iV(VS£ f /|| W) <: ll/H W) , 



and for the dual exponent p', 

(5-6) \\N,((sZW\W m ^ \\f\\iy m , 

(5-7) \\mVSZV)f\\iS £\\f\\is m - 

Analogous bounds hold for L*. The implicit constants in all inequalities and s > 
depend upon the standard constants only. 

Here and below 

(5.8) (S t Dj) f{x) := / -H-Elx, t; y, 0) f(y) dy, 1 < j < n, 

Jr" dyj 

(5.9) (S t D n+l )f(x)-- [ —E(x,t;y,s)\ s _ f(y)dy, 

JR" OS 
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and we set 
(5.10) 

n+l 

(S,V) := ( (JS,D{) , (S,D 2 ) , (S t D n+l ) ) , and (S f V) • /:= (S,D 7 ) / ; - , 

where / takes values in C" +1 . Similarly, 
(5.11) 

n 

(5 f V,|) := ((S^DO,^^),...,^,,)) , and (S t V {l ) ■/:=£) //, 

7=1 

where / takes values in C". 

Note that for ?-independent operators, we have by translation invariance in t that 
(S t D n+1 ) = -d,S,. 

Remark 5.12. We remark that, considering only the tangential gradient Vy in (5.6), 
the latter may be re-formulated as 

(5.13) sz ■ ill -» T£(K + l y 

We further note that (5.5), (5.6), ?-independence of L, and the De Giorgi-Nash- 
Moser bounds immediately imply that 

(5.14) 8 t Sf : L p (R n ) -> r£(Rf '), 1 < p < oo, 
under the assumptions of Proposition 5.1. 

Proof of Proposition 5.1. The square function bound (5.2) was proved in [R] (for an 
alternative proof, see [GH]). The fact that (5.2) implies (5.3) and (5.4) is basically 
a combination of results in [AAAHK] and [AA]. See Proposition 1.19 in [HMM] 
for a detailed discussion and references. The fact that (5.3) for L and L* implies 
(5.5) and (5.6) has been proved in [HMiMo, Theorem 1.1 and estimate (4.45)]. 
Finally, (5.7) can be found in [HMM]. □ 

Proposition 5.15. Suppose that L is t-independent, that L and its adjoint L* satisfy 
the De Giorgi-Nash-Moser bound (1.6). Then 

(5.16) VZT 1 div : T{ -» ?£, 1 < p < 2 + e, 

(5.17) VL- y \ :T[ ^fL, \<p<2 + e, 

for some s > 0, depending on the standard constants only. Here the operator L~ l y 
is to be interpreted via 

(5.18) (L' 1 !^) (y,s) := II Eiy^-xj^ixj)—, (y,s)eR n + +i . 

Proof. The Proposition was proved in [HMM]. □ 

Proposition 5.19. Suppose that L is t-independent, that L and its adjoint L* satisfy 
the De Giorgi-Nash-Moser bound (1.6). Then there exists e > such that for all 
1 < p < 2 + e, we have 

(5.20) VD L ±t : Li 7£(Rf '). 
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In particular, 

(5-21) l|tf((VD*/)|| W) < C||V||/|| W) . 

//ere C > ant/ e > depend upon the standard constants only. Analogous bounds 
hold for L. 

Proof. By the usual density considerations, it is enough to demonstrate (5.21) for 
/ e Cq(R"). To start, let us concentrate on the estimate for d t Df, that is, 

(5-22) \\N(d t Dff)\\ LP{Rn) < HVii/lboR*). 

To this end, we note that L* y s E*(y, 0; x, t) = for all (x, t) e R'_[ +1 , since (y, 0) is 
on the boundary. Hence, one can formally write 



(5.23) 



d t dy A ,, y E*(y,0;x,f)f(y)dy= d s e„ +1 A*(y)V y<s E*(y,s;x,t) 



f(y)dy 



J(y)dy 

s=0 



= - f J2 djejA*(y) V y ,, E*(y, s; x, t) 
- f L*E*(y,0;x,t)f(y)dy- [ Y,d J A* Jn1rl <y)d l E*(y,s;x,t) 

JR" JR" ~1 



f(y)dy 

s=0 



A;(yW y E*(y,0;x,t)V y f(y)dy 



r n 

+ / y^A* n+l iy)d s E*{y,s;x,t) djf(y)dy. 

JR" ~[ S=Q 



Of course, the formal computations in (5.23) should be interpreted in the weak 
sense, using, in particular, the weak definition of the normal derivative (2.10) (see 
the discussion in proof of Theorem 5.35 for more details). 

The desired estimate on the second term on the right-hand side of (5.23) follows 
from (5.14). Passing to the first term, recall that the tent space r£ can be realized 
as a space of the linear functionals on T p x , where the latter is defined as a collection 
ofFeM(Rf 1 )such that 



where 



WU ■= \\Ci(W 2 F)\\^ < oo, 
1 1 

W 2 F(x, t) = ( ff \F(y, s)\ 2 dyds J , (x, t) e Rf 1 . 

\W(x,t) I 

Indeed, it was proved in [HR], Theorem 3.2, that r£, = \T\ \ , 1 < p < oo. 
Pairing the first term on the right-hand side of (5.23) with O, one obtains 
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(5.24) ff ( Al(y)V y E*(y,0;x,t)V y f(y)dy<I>(x,t) 

JjR"+ l JR" 



dxdt 



^A;(y)V y T®(y,0)V y f(y)dy, 

where T - (L*)~ l ~, that is, as before, 

C* (~* d dt 

(5.25) T4>(y,s) = T s <I>(y):= E*(y,s;x,t)0(x,t) —. (y,s)eR n + +1 . 

JjR'l +1 t 

The goal is to show that 
(5.26) 



A;(y)V y TO(y,0)V y f(y)dy 



< WW?* HV||/|b>, 



for any O e T[ . 

It is sufficient to verify (5.26) for <I> smooth and compactly supported in R" +1 , 
as such functions are dense in Tf . At this point, recall the estimate (4.11) or, more 
precisely, (4.12). One can carefully track the proof of (4.11)-(4.12) to see that all 
the computations are justified for v := T§> with O smooth and compactly supported 
in R" +1 . However, v = TO is not a solution in R" +1 and hence, the term Illy (that 
is, the last integral on the right-hand side of (4.12)) will not be annulated. Instead, 
we have L*v(x, t) = L*T®(x, t) = 0(x, t)/t for (x, t) e R* +1 . 

All in all, it is enough to bound the right-hand side of (4.12) with v = TO, 
<D € ff(R" + +l ). Thus, one has to prove: 

\\M(t VdtTMWm*) * ^kfiRT'y I^*(W»IIw(r«) * PHff (K n y 

or, equivalently, in the language of tent spaces, 

(5.27) \\mTM T(m+1) < Pll^i,. 

(5.28) IWr^H^i, * Pllsf^i,. 

and to bound the last term on the right-hand side of (4.12) with / in place of /. 
However, 

(5.29) II d 2 P t f(x) L*T<b(x,t) tdtdx = II d 2 t PJ(x)<5>(x,t)dtdx, 
J J • JJ : 1 

so that using once again duality relationship for tent spaces the desired bound on 

(5.29) reduces to 

(5.30) \\N(tdjP t f)\\ Lmn) < IIVy/lbcR*) = ||V||/|| W) . 

The estimate (5.30) follows from the bound N(td 2 P t f) £ M(V\\f), while the latter 
can be proved essentially by the same argument as that in (4.6)-(4.9), using the 
off-diagonal decay estimates on t 2 d 2 P t and the Poincare inequality. 

It remains to discuss (5.27) and (5.28). A direct computation shows that the ad- 
joint of the operator d t T, under the usual tent space pairing (O, f) = Jj^n+i O f 

is -d s L~ l j (since by ^-independence one can swap the derivatives in t and s on the 
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fundamental solution). Since T^, = (rf \ , 1 < p < oo, (see, e.g., [HR], Theorem 
3.2, or the duality argument in [CMS]), (5.28) follows from (5.17) in the range 
1 < p < 2 + £. 

Finally, the adjoint of the operator tVd t T, under the tent space pairing (<!>, *F) = 
ff Rn+l ^f-, is -d s L~ l div, and T{ = fe'V, 1 < p < oo, (see [CMS]), so that 
(5.27) follows from (5.16), once again, in the range 1 < p < 2 + e, as desired. 

This finishes the proof of (5.22). 

We claim that the desired bound (5.21) follows from (5.22), or, to be precise, 
from the estimate 

(5.31) WNMDfDWwm.n) < QWufW^ny 

Since d t D\f is a solution, (5.31) is equivalent to (5.22) by the De Giorgi-Nash- 
Moser estimates. 

Indeed, using the argument in [KP], p. 494, we see that the left-hand side of 
(5.21) is controlled by the sum of ||M(V||Df 
where M is the Hardy-Littlewood maxima 
terms is bounded by ||V||Z)f 
side of (5.31). 

Thus, it remains to show that 



i=o /)IIzar») and \\M(N*(d t Dtf))\\iJ>(M»)> 
function on R". The first of these two 
f\\ip- The second one is bounded by the left-hand 



(5.32) V,|Df 



: L?(R") Z/(R"), 1 < p < 2 + s. 



Let /, g € C^°(R"), (where g is vector- valued) and set w(x, t) := St* t divy g(x), 
x e R", so that, by (5.6), 



(5.33) M-MWmn^WgWv 



Note that w is a solution (for L*) in the lower half-space, and choose an appropriate 
solution v in the lower half-space such that d t v = w. Dualizing, and using the 
equation, we then have 

(5.34) / V||D f f{x)g{x)dx= / f(x)d VA ,w{ x ) dx = / /(*) d t d VAm v(x) dx 
= -Y,f djf(x)A* n+l d t v(x)dx- [ V||/A*Vdx =:/ + //. 

~[ JR" JR" 

For term /, we just use (5.33), and for term //, we apply estimate (3.9) and then 
(5.6) and (5.7) (one can go over the argument of (3.9) and justify integration by 
parts and convergence of involved integrals even though the data of w on the bound- 
ary is not Cq). □ 

Theorem 5.35. Let L be an elliptic operator with t-independent coefficients such 
that the solutions to Lu = and L*u = in R" +1 satisfy the De Giorgi-N ash-Mo ser 

estimates. Let u he the solution to the Dirichlet problem Lu = in R" + , u = f, 
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for some f e C™(RP), in the sense of Lemma 2.4. Then for 1 < p < 2 + e 

(5.36) l|W(Vw)||£P(Kn) < C (||V||/|| W) + \\d VA u\\ mK r } ) , 
with C and e depending on the standard constants only. 

Proof. By Green's formula, the solution to Lu = in R" +1 , u = f, can be 

J J + dR1 +[ 

written as 

(5.37) u(x,t) = - d VA ,, y E*(y,0;x,t)f(y)dy+ E(x,f,y,0)d VA u(y)dy, 

JR n JR" 

or, in the language of layer potentials, 

(5.38) u = -Dt(f)+Sf(d VA u). 

With (5.38) at hand, (5.36) follows directly from (5.21) and (5.5). 

It remains to justify the use of the Green's formula in the present context. 

To this end, fix X := (x ,t ) e R^ +1 , and let e > 0, with s < t /S. We 
define a nice approximate identity on in the standard way as follows. Let 
CD £ C^(R" +1 ), with CD > 0, supp(cD) c 5(0, 1), and JJ, n+] G> = 1, and set cD e (X) := 
£-"- [ ®(X/£). For (y, s) e R" +1 , we set 



H £ (y, s) := (O e * E*(y, s, ■, -)j(x , h) . 
We observe that by definition of the fundamental solution E* , 

where Of°(X) := O £ (X - X) e C^(B(X , e)). In particular, then, we have that 
H E e W U2 (R n+1 ). 

Let F be a C ( ~(R' ,+1 ) extension of / = u(; 0), and set 

u(x, t), t>0 
F(x,t), t<0. 



u(x, t) := 



Since e <sc to, and since L*H E = Of in the weak sense in R" +1 (and in particular, 
L*H £ = in the lower half-space), we have 

(5.39) ® £ *u(X )= \ \ u(X)®x°(X) dX - ff Vu ■ A*VH £ 

= ff AVuVH~ £ +ff VF ■A^VJT, - (d VA u,T £ ) - (b\jT £ , f), 

JJR" + +1 JjR n _ +1 

where h £ := H £ (y, 0), and where in the last step we have used Lemma 2.11. By 
definition of H £ , (5.38) follows, letting e — > 0, once we establish the convergence 
in e for both terms on the right-hand side of (5.39). 

The convergence of the first term is quite easy, as h £ — > E* in Ljy 2 (R"). 

Let us discuss the second one. In this regard, one has to show that for e <sc 
to, the variational co-normal derivative -d VA ,H £ , initially defined in the sense of 
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Lemma 2. 1 1 , belongs to L 2 (K"), and satisfies 

(5.40) -d VA .H E (y,0) = e n+1 ■ A* (® e * V y , s E*(y, s, ■, •) | s=Q j (Xq). 

That the right hand side of (5.40) belongs to L 2 (R n ) follows directly from [AAAHK, 
Lemma 2.8]. Thus we need only verify the identity (5.40). By Lemma 2.11, ap- 
plied in the lower half-space, it is enough to verify that for any F e C ( ~(K.' l+1 ), with 
F(-, 0) := /, we have 

(5.41) fJ^VFA*W s = j J(y)rf -A*(® e *V y<s E*(y,s,-,-)\ s=0 )(Xa)dy. 

To this end, let P n be a nice approximate identity in R", applied in the y variable. 
By the divergence theorem 



/ /O0 $ ■ P v U* (o £ * V y , s E*(y, s, ■, -)\ ) (Xo)) dy 

JR" V V ■ / / 

-ff 



div F(y, s) P„ (A* O e * V y , s E*(y, s,;-)) (X ) dyds 



ff 



VF(y, s)-P„(A*VH £ (y, s))dyds, 



where in the last step we have used the definition of H e , and that it is a null solution 
of L* in R n _ +1 . We now obtain (5.41) by letting 77 -» 0. 

This finishes the argument. □ 

Corollary 5.42. Assume that L is an elliptic operator with t-independent coeffi- 
cients, and that L and L* satisfy the De Giorgi-Nash-Moser condition. If for some 
1 < p < 2 + e the Dirichlet problem (D p r)for L* is solvable in R" +1 and, in ad- 
dition to (1.4), the solution satisfies the square function estimates (1.12), then the 
Regularity problem (R p )for L is solvable in W£ l . That is, assertion (a) implies (b) 
in Theorem 1.11. 

Proof. The combination of Theorems 5.35 and 3.1 entails that every solution to 
Lu = in R'i +1 , u = f e C^°(R re ) in the weak sense of Lemma 2.4 satisfies the 

estimate (1.5). With this at hand, the solvability for a general elliptic operator L, as 
defined before the statement of Theorem 1.11, follows immediately. □ 

6. Invertibility of layer potentials. Layer potential representations of 

solutions 

Propositions 5.1 and Remark 5.12 ascertain that the layer potentials are always 
bounded for elliptic operators with ^-independent coefficients that satisfy the De 
Giorgi-Nash-Moser bounds. We note that, furthermore, (5.5) and Lemma 2.2 imply 
that for / e Z/(R"), 1 < p < 2 + e, the single layer potential S L f has well-defined 
normal and tangential derivatives on the boundary in the sense of Lemma 2.2 (de- 
noted below by -e n -A VS L f and V||5 L /, respectively). Both e n -A VS L f and VyS/ 
belong to LP with the appropriate estimates. In particular, 

(6.1) S L : L p (R n ) -> Lf(R n ), 
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is a bounded operator for all 1 < p < 2 + e. When necessary, we shall use the 
superscripts + and - to underline the limit taken from the upper and the lower half 
space, respectively. 

Recall the following jump relations: 



The jump relations (6.2), (6.3) can be found, e.g., in the proof of [AAAHK, 
Lemma 4. 18] for g € L 2 (R"), and then extended to any g € L P (MP), 1 < p < 2 + e, 
using (5.5) and Lemma 2.2. In particular (S L ) + = (S L )~ regarded as operators 
L P (W) -> L\ (R"), 1 < p < 2 + s. Similarly, (6.4) is established for / € Cg°(R n ) in 
[AAAHK, Lemma 4.18] and can be extended by continuity to / e Lj(R") for all p 
such that (5.21) holds. 

Proposition 6.5. Assume that L is an elliptic operator with I -independent co- 
efficients, and that L and L* satisfy the De Giorgi-Nash-Moser condition. As- 
sume further that for some 1 < p' < oo the Dirichlet problem (D p >) for the 
operator L* is solvable with the square function bounds, that is, the solution u 
satisfies both (1.4) and the square function estimate (1.12). Then the operator 
(5' i ) ± : L P (R") — > L p (R n ) is invertible, - + y = 1. An analogous statement holds 



Proof. First of, the boundedness of (S L ) ± follows from the discussion of (6.1) 
above. The bounds from below can be established as follows. By (6.2), (6.3) 
and Lemma 2.2, for any g e C^(R") we have 

(6.6) \\ 8 \y m = \\i n ■ A<ys L ) + 8 - i„ ■ A^s L yg\\ mv) < c ii(v B s Vibw. 

To justify this, we point out that whenever g e C^°(R n ), the functions (S L^g belong 
to the space Lj(R n ) n Z/(R") n L«(R n ) from Lemma 2.4 and to t\ (W). This can be 
seen from (6.1) and its dual. Hence, (S^g can be approximated by fy e C^°(R") 
in Lj(R") n L r (R n ) n L«(R") and in t\ (W). Using Lemma 2.5 in [AAAHK], one 
can also show that the solution S t g belongs to W 1,2 (R" +1 ). Now, for each such 
fk there exists a weak solution of the Dirichlet problem by Lemma 2.4 (call it u^) 
and, due to (2.5), these solutions converge to S t g in W 1,2 (R") norm. Moreover, Uk 
satisfy Theorem 3.1. Hence, the corresponding weak conormal derivatives d VA Uk 
form a Cauchy sequence in LP and thus, converge to some L p function, h, in L p 
norm. On the other hand, as mentioned above, converge to S t g in W 1 ' 2 (M n ) 
norm and hence, d VA Uk converge to the conormal derivative of S t g in the sense of 
distributions. Thus, the conormal derivative of S t g can be identified with h and its 
L p norm is bounded by the LP norm of the tangential derivative, as desired. 

Next, we note that (6.6) for all g € C^°(R") implies (6.6) for all g e LP(R n ). In- 
deed, every g € Z/(R") can be approximated in L p by functions g%, and by (6.1) 
(V\\S L g) ± is well-defined, belongs to L p , and the sequence (V\\S L gk) ± converges in 
LP to (Vi^^)*. Thus, (6.6) holds for all g g L P (R") as well. 



(6.2) 
(6.3) 
(6.4) 



i n ■ A(vs L y g - i n ■ A{vs L y g = -g, 
(v ll s L ) + g-(v ll s L y g = o, 
(D L rf - (D L rf = -f. 



forV. 



30 



STEVE HOFMANN, CARLOS KENIG, SVITLANA MAYBORODA, AND JILL PIPHER 



Let us show that (S 1 )* is surjective 2 . Choose some / € C™(W). Let vr be the 
solutions to the Dirichlet problem in the weak sense of Lemma 2.4 with boundary 

data / in R" +1 . By Lemma 2.2, the distributions g ± = e n ■ AVw* 
We shall show that 



lie in WR"). 



(6.7) (S L ) ± (g.-g + )=f in Lj(R"). 

To this end, observe first that are also the solutions to the regularity problem 
(R p ) (see, e.g., the proof of Corollary 5.42) and moreover, according to the proof 
of Theorem 5.35, we have 

u ± = -(D^) ± (f) + (S^) ± (g ± ). 
Using jump relations (6.2)-(6.4), we deduce that 



- u 



f+{s L n g+ -g„). 



Restriction to the boundary is interpreted here in the sense of the non tangential 
limit of the gradient, and the equality holds in L^(R n ). We used, in particular, the 
fact that (S L ) + - (S L )~ regarded as operators LT(R n ) -> L\ (R"), alluded to above. 
On the other hand, by definition of u as a solution to (R p ) with data /, we have 

— u~ = 0, again in the sense of the non-tangential limit of the gradient and 

R" 

equality in L\ (R n ). The combination of these two facts immediately yields (6.7) 
for / e C^(R n ). Now, for general / € L\ (R") we can use the limiting procedure, 
(6.7) for functions and (6.6) to find g e L P (R") such that (S^g = f. □ 

Corollary 6.8. Assume that L is an elliptic operator with t-independent coeffi- 
cients, and that L and L* satisfy the De Giorgi-Nash-Moser condition. 

If (R p ) is well-posed in K.+ +1 for some 1 < p < 2 + e then the solution can 
be represented by means of layer potentials, that is, assertion (b) implies (c) in 
Theorem 1.11. An analogous statement holds for L*. 

Proof. The corollary is a combination of the boundedness results in Proposition 5.1 
and invertibility of layer potentials following from Proposition 6.5. □ 

Corollary 6.9. Assume that L is an elliptic operator with t-independent coeffi- 
cients, and that L and L* satisfy the De Giorgi-Nash-Moser condition. 

If (Rp) is well-posed in R+ +1 for Lfor some 1 < p < 2 + s with the solution 
represented by means of layer potentials then (D p >) is well-posed in R" +1 for L*, 
with the solution represented by means of layer potentials, that is, assertion (c) 
implies (d) in Theorem 1.11. An analogous statement holds for L*. 

Proof. By Proposition 6.5 and preceding discussion, the operator (5' L ) ± : L p (B, n ) — > 
Lj (R") is bounded and invertible. It follows that the operator 

< L*\± . T p' 



(S L ) ± : L^(R") -> V (R n ) 



We note that the surjectivity argument presented here was inspired by the surjectivity result for 
layer potentials in [BM] 
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is bounded and invertible. Taking in account (5.6), the solvability of the Dirichlet 
problem via the representation (1.15) follows immediately once we observe that 
h e Z/j if and only if there exists L4 e LP' such that divy fl = h, with the accompa- 
nying norm equivalence. □ 



7. Well-posedness and uniqueness 

Sections 3-6 establish equivalence of solvability of the Regularity and Dirichlet 
boundary value problems, in the context of Theorem 1.11. However, even having 
solvability at hand, well-posedness turns out to be a delicate issue. Basically, under 

1 2 

current assumptions on operator, the notion of a solution as a W^ c function satis- 
fying appropriate non-tangential maximal function estimates does not seem to be 
sufficient to justify, e.g., the Green's formula. Thus, one has to appeal to a slightly 
weaker notion of uniqueness and, respectively, well-posedness of (D p >), (R p ) in the 
current context. 

Let us first note that for every / £ Cq(R") the solutions to (D p >), (R p ) obtained in 
Theorem 1.11 are the weak solutions in the sense of Lemma 2.4. And, in the class 
of weak solutions, they are unique (again, by Lemma 2.4). Clearly, it is sufficient 
to justify this for layer potential representations (1.14), (1.15). To this end, observe 
that for / e Cg°(R") we have 

VSt(Sfr l fe f| TP, 

l<p<2+£ 

~ In 

byTheorem 1.11. Thus, in particular, VS^Sq) -1 / £ T^ +[ and hence, by Lemma 2.2 
in [HMiMo], we have V5f(5^) -1 / £ L 2 (R" +1 ). Similarly, 

S^S^fe f| Tt => SHsfr l fe f| L^Rf 1 ) 

l<p<2+£ \<p<2+E 

implies that the 1/(1 + |X| 2 )-weighted L? norm of SfiS^)' 1 / is finite. Altogether 
then, SfiSft)- 1 / e W 1 - 2 ^ 1 )- 

Having this information at hand, we can define the solution to (D p >) (resp., (R p )) 
with arbitrary data / £ LP (resp., / € L p ) as a limit. That is, one takes for any / £ 
LP' a sequence {fk)k of Cq functions approximating / in LP' norm. Corresponding 
to each one can find a unique weak solution £ W l ' 2 (Wl +[ ) with the data 
on the boundary. Then, if (D p >) is solvable, form a Cauchy sequence in T£,. 

The limit of this sequence in TL is the solution u with data /. By definition, u 
satisfies the desired estimates, and the standard argument shows that two different 
sequences converging to / yield the same u. One moreover can easily verify that 
the solution obtained as such a limit also satisfies the layer potential representation 
(1.15) as Sf(S^r l f k converge to SfiS^)' 1 / in the norm of t£. 

Recall that the solutions in Theorem 1.11 were already defined by this limiting 
procedure (see the Remark after Lemma 2.2). The difference is that now we specify 
Uk for every in a unique fashion. And hence, the solution u defined as a limit of 
such Uk's is unique and the problem is well-posed. The same argument applies to 
the Regularity problem. All in all, both in the Dirichlet and Regularity case, one 
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has uniqueness of solutions (in the class of the limits of the weak solutions in the 
appropriate tent space) and thus, one obtains well-posedness in this sense. 

In fact, we shall prove a stronger result. Recall that (D p >) (resp., (R p )) is well- 
posed with the Green formula representation if for every / € LP (resp., / e Lj) 
there exists a unique solution u satisfying (1.4) (resp., (1.5)) such that the Green 
formula (5.37) holds for u. Let us single out the following result. 

Lemma 7.1. Under the conditions of Theorem 1.11 the corresponding boundary 
problems are well-posed with the Green formula representations. In particular, 
uniqueness holds in the following sense. Take the operator L falling under the 
scope of Theorem 1.11 and, respectively, 1 < p < oo such that the statements 
(a) — (d) are valid. Then whenever u € VK/ '^(K." +1 ) is such that N*u € LP (K.' ! ), 

and the Green formula (5.37) holds, we have u = 0. Similarly, whenever 



an 



u € Wj£(R n + +] ) is such that AT(Vh) e Z/(R"), u 
(5.37) holds, we have u - 3 . 



= and the Green formula 

da 



Proof. Indeed, first of all, as we discussed above, for every data the solutions 
provided by Theorem 1.11 are the weak solutions. The proof of Theorem 5.35 
demonstrated, in particular, that for every weak solution the Green formula rep- 
resentation (5.37)-(5.38) holds. Using mapping properties of the layer potentials 
(5.21), (5.5), and (2.3), one can show now that any solution to (R p ) defined by the 
limiting procedure as above also satisfies (5.37)-(5.38). 

Similarly, using (5.6), we see that any solution to (D p >) defined by the limiting 
procedure as above satisfies (5.37)-(5.38) provided that d VA u k converge to d VA u in 
L^(R n ). However, u k = Sf(S^y l f k for all k e N and, as we discussed above, u 
defined by the limiting procedure necessarily satisfies u = Sf(So) -1 /- Thus, one 
just has to show that d VA S^ is a bounded operator in L p _ x . This, however, is 
exactly the dual statement to (5.32). 

It remains to show that solutions satisfying the Green formula (5.37)-(5.38) are 
unique. To this end, assume, for instance, that u solves (R p ) with / = and satis- 
fies (5.38). Then u = S^(d VA u). Taking the non-tangential limit on the boundary 
(by Lemma 2.2) and using invertibilty of the single layer potential established in 
Proposition 6.5, we deduce that d VA u = 0. Hence, u given by (5.37)-(5.38) is 
identically zero, as desired. 

Turning to the Dirichlet problem, we first have to explain what is meant by a 
solution with the Green formula representation, or more precisely, how the normal 
derivative of an arbitrary u e W^(K." +1 ) with Lu = and N t u e L P '(W) is defined. 

Choose 1 < p' < oo such that {D p >) is solvable with the square function bounds, 
and assume that u e W^(R" +1 ) with Lu = and N*u € Z/(R") is a solution with 



data u 



= f£LP. 



3 The data of the Regularity problem are defined modulo constants (as elements of L^), and as 
usual, the uniqueness is understood in the sense of working with a particular realization of the data 
equal to zero a.e. 
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ip in the sense of 



Let us define a tempered distribution d VA u as 
(7.2) (d VA u,<p) := {f,dy A .<!>), V<pe q°(R"), 

where cp is the unique weak solution to L*cD = 0, ® 
Lemma 2.4. Note that the definition coincides with the one used before when, 

e.g., u is the unique weak solution to Lu = 0, u = f in the sense of Lemma 2.4. 

R" 

Let us show that d VA u, defined by (7.2), belongs to L^(R n ) = (Lj(R"))*. Indeed, 
we have taken 1 < p' < oo such that (D p >) is solvable. Thus, by Theorem 1.11, 
the Regularity problem (R p ) is solvable, and moreover, as we discussed above, for 
every ip e Cq the solution is given by the unique weak solution in the sense of 
Lemma 2.4 and thus, cp in (7.2) satisfies the estimate 

II<9 Va ,,®||lp(r<') ^ I Ml^gr) > 
by Theorem 3.1. It follows that 

\(d VA u,(p)\ < QlfW^^WipWiP^ny 

and hence, d VA u € Z/j(R"). The remainder of the uniqueness argument follows the 
same lines as that for the Regularity problem. □ 

Remark. Concluding this Section, we note that whenever L and L* satisfy the 
De Giorgi-Nash-Moser bounds in R" and at the boundary, the Regularity and 
the Dirichlet problems are well-posed in the usual strong sense (without explicitly 
requiring validity of the Green formula). That is, whenever u e w/ ' c 2 (R" +1 ) is 

- 0, we have u = 0. And similarly, whenever 

dn 

= 0, we have u = 0. This is the 

da 



such that N t u € LP (R re ), u 



u e W}^(R n + +l ) is such that 7Y(Vw) £ L P (W), u 
case, e.g., when the elliptic operator has real coefficients. 



8. Proofs of Corollaries 1.16-1.19 

We conclude the paper with the Proofs of Corollaries 1.16-1.19. 

Proof of Corollary 1.16. The De Giorgi-Nash-Moser bounds are stable under L°° 
perturbation of the coefficients. Thus, if Lq falls under the scope of Theorem 1.11, 
so that Lo and L ( * } both have the De Giorgi-Nash-Moser property, then so do L and 
L*. Thus, the conditions (a) - (d) of Theorem 1.11 are equivalent for the operator 
L too, and it is sufficient to prove one of them. It is easier to access (c) or (d). Let 
us focus on (c). 

The fact that / ^ N(VS^f) is bounded in Z/(R"), 1 < p < 2 + e, for any 
£-independent operator L (with some e depending on the ellipticity constant) is 
known - see the discussion of (5.5). Respectively, 

: L p (R n ) ^ Lf(R M ), 

is bounded by Lemma 2.2. The operator norm in both cases depends on standard 
constants only. 
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As far as invertibility is concerned, we have by analytic perturbation theory (see 
the Appendix for details) 

(8-1) II^II^O ~~ ^11^0 \\LP(R n )^LP(R") ^ C||A - A|| L o=( R n), 

and hence, invertibility of S q follows from that of S q° via the Neumann series. □ 

Proof of Corollary 1.17. This is essentially the main result of [HMM]. Let us 
provide some explanations. 

It is clear that it is sufficient to prove (c) and (d) for L, together with the square 
function bound (1.12) for the solution given by (1.15). Then (a) and (b) would 
automatically follow. We note that equivalence of (a) - (d) for L in this case does 
not follow from Theorem 1 . 1 1 , as the coefficients of L are not ?-independent. 

The fact (c) and (d) for Lo imply (c) and (d) for L under the conditions of the 
Corollary is one the main theorems in [HMM]. □ 

Proof of Corollary 1.19. Recall that the De Giorgi-Nash-Moser bounds for solu- 
tions inside the domain and at the boundary are always valid for the operators with 
real coefficients. 

The validity of (a) for an elliptic operator with real, f-independent, possibly 
non-symmetric, coefficients is the main result of [HKMP]. Then the validity of 
all ascertions (b) - (d) of Theorem 1.11 for such an operator follows from The- 
orem 1.11, and the claimed perturbation results follow from Corollaries 1.16 and 
1.17. □ 



9. Appendix: Analyticity of VS t 

In this section we present a discussion of (8.1). It is essentially treated within the 
framework of the analytic perturbation theory in [Ka] (and was already used, e.g., 
in [AAAHK], [B]). However, a detailed argument for the particular case at hand 
does not seem to be available in the literature, and for completeness, we provide it 
here. 

Let Ao be elliptic, (n + 1) X (n + 1), ^-independent, and bounded measurable, 
and let A z := Aq + zM, where M is (n + 1) X (n + 1), ^-independent, and bounded 
measurable, with ||.M||z«>(r») < 1- Set Lo := -divAoV, and L z :- -divA z V, and 
suppose that null solutions of Lo satisfy the De Giorgi-Nash-Moser bounds. For 
|z| small enough, say |z| < sq, we have that L z is also elliptic, and satisfies the De 
Giorgi-Nash-Moser bounds. By ellipticity, for |z| < £o> 

VLJ 1 div : L 2 (R" +1 ) -> L 2 (R' ,+1 ), 

or equivalently, 

L7 1 : W~ h2 (R n+1 ) W l ' 2 (R n+l ) . 
Moreover, we have the Taylor expansion 

DO 

VLJ 1 div - VLq 1 div ^ (zMVLq 1 div)* , 

k=Q 

which is convergent, as a mapping from W~ l ' 2 (W +l ) to W l ' 2 (W +l ), if eo is small 
enough. Therefore, the mapping z — > L" 1 , taking values in the space of bounded 
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operators from W~ l ' 2 (M. n+l ) to W l ' 2 (M. n+l ), is analytic in a neighborhood of z - 0. 
The same is true for L*, so by trace theory, we then have that z — > Tro(L*) _1 
is an analytic mapping, taking values in the space of bounded operators from 
W~ l ' 2 (R n+1 ) to H 1/2 (R n ). Here Tr denotes the trace operator, on R"x{0} - d(R'j +1 ). 
We then define the single layer potential for L z , denoted by S Lz , as the adjoint of 
the operator Tr o(L*) _1 , so that z — > S L '~ is an analytic mapping taking values in the 
space of bounded operators from H~ l ^ 2 (W) to W 1,2 (K" +1 ). By trace theory again, 
we have that z — > S f : is an analytic mapping taking values in the space of bounded 
operators from H~ l ^ 2 (M. n ) to /f 1,2 (R"), where sf z denotes the restriction of S Lz to 
the hyperplane x n+ \ = t; i.e., z — > VyS; is an analytic mapping taking values in 
the space of bounded operators from //~^ 2 (]R") to H~ l ' 2 (W). Thus, 

z^<V||5^/,g) 

is an analytic function, for every / e H~ l ^ 2 (R"), and every g e //^ 2 (]R"), in partic- 
ular, for every pair f,g€ C^(K"). By [HMiMo, Theorem 1.1], we have that 

sup sup ||V||5 , f z || L p(R«)^L"(R") <C P , 1 < p <2 + e, 

k\<£0 * 

with C p depending only on p, ellipticity, dimension, and the the De Giorgi-Nash- 
Moser constants. Therefore by [Ka, p. 365], since C$ is dense in LP and LP , we 

have that for each fixed t, z — > ^\\Sf z is an analytic mapping, taking values in the 
space of bounded operators on L p (E, n ), in the disk |z| < so. This means that 

II d t I, 

so that 

(9.1) sup ||V||5f ; - V||Sf°|| L , w ^ L , w < C p |z| . 

t 

Now, given A as above, takeM :- (Ai-A )/||Ai-A ||l"(R")> so that A\ = A +z\M, 
with z\ = |Ai - A ||l~(r»)- If l|Ai - A ||l»(r») < £o, by (9.1), we have that 

SUp llVySf 1 - V||Sf°||LP(R»)-»LP(R") ^ C p ||Ai - AoI|l°°(R") , 

as desired. 
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